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Abstract 



O 
O 

o 

<. 

, Consequences of different discretizations of the two-dimensional Dirac 

operator on low energy properties (e.g., the number of nodes) and their 
relations to gauge properties are discussed. Breaking of the gauge in- 
J> ■ variance was suggested in a recent work by M. Bocquet, D. Serban, 

and M.R. Zirnbauer |]cond-mat/9910480| in order to destroy an inter- 
ff^ I mediate metallic phase of lattice Dirac fermions with random mass. It 

I is explained that such a procedure is inconsistent with the underlying 

■ lattice physics. Previous results point out that the logarithmic growth 

■ of the slope of the average density of states with the system size, ob- 
, tained in the field-theoretical calculation of M. Bocquet et al., could 

^1 be a precursor for the appearence of an intermediate metallic phase. 

o 
o 

Two-dimensional (2D) Dirac (Majorana) fermions can be derived in statistical 
^ ! physics (2D Ising model [|I|) and in models of condensed matter systems (quasi- 

^ [ particles in d-wave superconductors in the resonant valence bond state of the 

two-dimensional Heisenberg model 0], or in the quantum Hall effect ) al- 

ways from lattice models. The relation of all these lattice models and 2D Dirac 
fermions is based on the agreement of their low-energy properties on certain length 
scales. Using continuous Dirac fermions, which have a linear dispersion, their k = 
node is identified with one of the nodes of the original lattice model. However, 
there are technical reasons to study a discrete (lattice) version of the Dirac operator 
because of problems related to the unrestricted spectrum of the continuum model. 
In the presence of randomness there is also the particular problem of correlations 
on a characteristic length scale which restricts the use of the unrealistic white noise 
randomness. In the renormalization group approach or in the saddle-point ap- 
proximation it is sufficient to specify the lattice by a cutoff of the wavevector, 
assuming that only the k = node is important. For a more detailed discussion of 
the model one must define the lattice structure explicitly. The simplest case would 
be a nearest neighbor (NN) discretization of the lattice difference operator 



I. INTRODUCTION 
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V,/(r) = -[/(r + e,)-/(r-e,)] 



where Cj is the lattice unit vector in direction j on a square lattice. This choice, 
which was used in several papers by the author, exhibits an "intermediate phase" , 
characterized by a non- vanishing density of states (DOS), in the presence of a ran- 
dom mass for < nic- It turned out that this phase is metallic with finite 



conductivity [JI2 



The discretization ([l|) was criticized in a recent work by Bocquet et al. [0 



because of the associated sublattice symmetry. It was argued that the non-vanishing 
average DOS of this model is a sole consequence of the fact that it belongs to a 



special universality class, and that a different discretization would have a vanishing 
average DOS. Their argument is based on the renormalization group calculation of 
Ref. which gives a linearly vanishing DOS at = 0. (However, problems with 
this calculation exist because the RG has a strange behavior in which the slope of 
the DOS diverges logarithmically with the system size |jl3|-) Since the non- vanishing 
DOS represents an intermediate metallic phase, observed and discussed for several 



models in Refs. P JT^JT^J15|JT6|JT7|] , the change of the discretization would result in a 
destruction of this phase. 

In the following, an interpretation of the 2D Dirac fermions with discretization 
(111) is given in terms of a tight-binding model with flux vr. It will be shown that the 
additional sublattice symmetry is a special gauge symmetry. The destruction of the 
gauge symmetry by an additional next-nearest neighbor term in @) leads to a lattice 
model with an inconsistent gauge field representation. More general, a change of 
the discretization affects physical quantities like the Hall conductivity which depends 
on the number of nodes. Finally, it will be shown, using a result from a previous 
calculation, that breaking of the sublattice symmetry does not neccessarily result in 
a vanishing average DOS at ii^ = 0. This indicates that even though the properties 
of the intermediate metallic phase vary with a change of the discretization, its very 
existence is rather stable. 



II. DISCUSSION OF THE MODEL 

The 2D Dirac Hamiltonian is given as 

H = Ho + mas, Hq = iVicri + iV20"2 

with Pauli matrices aj and a random mass with (m^) = 0. Vj can be any discrete 
lattice difference operator. It must be antisymmetric in order to have a Hermitean 
matrix H. For instance, it could be a NN difference operator (|I]) or a combination 
of a NN and a NNN difference operator. The former will be called "tt flux dis- 
cretization" and the latter "flux breaking (FB) discretization" for reasons which are 
explained subsequently. The additional NNN difference operator adds new nodes to 
the dispersion of the pure system. To understand the difference of the discretizations 
one can consider the pure model, i.e. the Hamiltonian without the random mass. 
Then the dispersion for the vr flux discretization is 
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Ej, = ± Y sin^ ki + sin^ k2 
which has nodes at kj = 0, ivr. For the FB discretization the dispersion is 



Efb = ±\I [sin ki + sin(2A;i)]2 + [sin k2 + sin(2A;2)]2. 

with nodes at kj = 0, ±7r, arccos(— 1/2). Additional nodes of the pure sytem can 
increase the DOS of the disordered system. This gives rise to the question whether 
all nodes contribute to the DOS at low energy or only some. On the other hand, 
the nodes with non-zero k can be removed by adding a diagonal hopping term to 
the Hamiltonian. As it was shown in a previous study [jlT| (and will be discussed in 
Sect. Ill), the k = node is sufficient for a non-zero average DOS at E = 0. Thus 
the states with wavevector k ^ appear with non-vanishing density. 

Two sublattices can be constructed by considering two layers of the original 
square lattice, where the layers refer to the Dirac index a = 1,2. Then the sublattices 
are defined through the condition that for the site (r, a) either ri -|- r2 + a is odd 
(sublattice Ai) or even (sublattice A2). The Hamiltonian with tt fiux discretization 
acts on these two sublattices separately: 

H = PiHPi + P2HP2, (2) 

where Pj is the projector on A^. Adding next-nearest neighbor (NNN) terms to (|I|) 
results in the Hamiltonian Hps which couples both sublattices. 

The transformation H HSa^ can be applied to the Hamiltonian with 



H commutes with So^ [|ri|,|I^. It changes the sign of P2HP2 but leaves PiHPi 
unchanged. This means it is a gauge transformation on A2. 

It is important to notice that property (@) does not mean that one can project 
on one of these sublattices and ignore the other projection because both sublattices 
are still statistically coupled through the same random mass. This can be seen in 
the product of the Green's function with its complex conjugate \{H + ^e)~J^,^, |^. It 



was shown [|T^ that for a' = a + 1 [mod 2) the symmetry of H implies 



and 

Thus I (if + i€)~ar'a'\^ '^^ '^^^ sublattice is the product of Green's functions from 
both sublattices. The non-trivial properties of the two-particle Green's function are 
a consequence of the Green's functions on the two sublattices which share a common 
random mass. 

The Hamiltonian Hj = PjHPj describes a tight-binding model on the sublattice 
Aj with hopping elements 
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Hrl,r±ei2 — -f^r2,r±eil — =t^ 
Hrl,r±e22 = -f^r2,r=Fe2l ~ =^1 

which represents a Hamiltonian with flux vr. Adding NNN terms to Vj changes 
this flux because the new terms have the same phase factors as the NN terms but 
on links which are twice as large as those of the NN hopping terms. This leads 
to an inconsistency in the model because particles would experience a magnetic 
field whose strength depends on the length of the hops. Moreover, the invariance 
property [H,S(r-^]^ = 0, which is the gauge invariance on sublattice A2, is violated 
by the NNN hopping term. 

It is possible to add NNN hopping terms to H which are consistent with the 
node structure of the vr flux discretization. The simplest is 

It preserves the gauge field of the NN term. Another one is 

which does not preserve the gauge field because it contributes a positive hopping 
element for a = 1 but a negative hopping element for a = 2. Nevertheless, both 
NNN terms obey the condition 

[hnnn, Sas]^ = Sa-s]^ = 0, (3) 

but only satisfies the defining symmetry of class D of Altland and Zirnbauer 

m 

11 h' ^ 

'''nnn ~ ~^l''nnn "^l- 

Property (H) can be used to estimate the non-vanishing average DOS at = 0, as 



demonstrated in Refs. 10,11 



The discussion of the DOS requires only the diagonal part of the one-particle Green's 
function Grr,aa = {H ± iecro)~r^aa5 which can be represented by a sum over closed 
random walks [T^ beginning at r and returning to it with the same Dirac index. For- 
mally, the random walk representation can be obtained from the hopping expansion 
of the Green's function 

{Ho + mas + ^eao)""^ = {mas + ieao)'^ [Ho{ma3 + ieao)'^ . 

l>0 

This expansion is convergent for sufficiently large e. After averaging over the random 
mass one can perform an analytic continuation to arbitrarily small e > 0. The loops 
experience an effective gauge field because of complex hopping elements. An example 
is a simple plaquette (Fig. 1). In general the flux per plaquette is = vr, in units of 
the flux quantum (pQ. An important property of the random walks follows from 

G{—m,ie) = a2G^{m,ie)a2 

because this means that the random walks are reversed by the change of the sign 
of the mass. An interpretation is that the currents in the model can be reversed by 
reversing the sign of the mass, which is related to the change of the sign of the Hall 
conductivity a^y with m 
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A. Hall Conductivity 



A quantity which is sensitive to the type of discretization but robust against disorder 
is the Hall conductivity cTxy It can be measured as the linear response to an external 
gauge field [^. In a pure system (m = const.) with Hamiltonian H = ma^ + hiai + 
h2<72 it reads [H 

m r 1 (Pk 

"""^ " 2" i (m2 + hl + hlfl^ (27r)2 

in units of /%. For m ~ only the nodes contribute significantly to the intergral: 

Tfi 

(Txy{m) = — /(|m|), 
\m\ 

where /(O) ^ is proportional to the number of nodes. In the case of the vr fiux 
Dirac operator this gives just a Hall step for each sublattice which means a single Hall 
step for the corresponding vr fiux tight-binding model. For the modified Hamiltonian 
Hq + ma-i + h'^nn Hall conductivity is 

1 r m + sin^ ki + sin^ k2 dPk 

2 J [{m + sin^ ki + sin^ /ca)^ + sin^ ki + sin^ fcs]^/^ {2%^ 

which gives 

1 TTl 

Cxy ~ const. H — - (m ~ 0). 

271 \m\ 

Thus, the additional term h'^^^ contributes a constant Hall conductivity at small 
values of m. These results indicate that approximations of the original model should 
not change the number of nodes nor violate the structure of the Hamiltonian by 
adding new terms of the type h'^nn order to get the correct Hall conductivity. 

A general rule is that all nodes of the lattice model should be included in the 
calculation of the low energy properties. This has severe consequence for the renor- 
maliztion group (RG) calculation which usually deals only with one large length 
scale. The various nodes (different length scales) will create additional couplings 
under the RG transformation, leading eventually to a strong coupling behavior. 



III. MODEL WITH SINGLE NODE 

To study the contribution of the k = node to the average DOS, all nodes of the vr 
flux discretization must be removed except for that at k = 0. This can be achieved 
by adding a diagonal term to H 

H^ = H+{A- 2)(T3 
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with2A/(r) = /(r+ei)+/(r— ei)+/(r+e2)+/(r— 62). Obviously, this Hamiltonian 
violates the sublattice gauge-invariance condition: 

Aa^S = —a^SA. 

Without randomness (i.e. m = 0) the dispersion is 



Ei{ki, ^2) = ±Y (cos ki + cos k2 — 2)^ + sin ki + sin k2 

which is non-zero except for the node at kj = 0. Now this Hamiltonian can be 
coupled by a random field m' to the Hamiltonian with nodes at ±7r 

H2 = H-{A + 2)a3 



as 



H 



Hi m'cx^i 



Both Hamiltonians Hi and H2 violate the sublattice gauge-invariance condition 
but preserve the symmetry of class D 



H = -aiH^ai. 



The Green's function can be projected on the subspace of Hi 

[Hi m'a^{H2 + ie)"V3m']~^ 

It was shown that this expression leads to a non-vanishing average DOS at = 



11], indicating that the projection on one node produces already a non-vanishing 



average DOS. 

As a remark it should be noted that the pure 2D Ising model at the critical point is 
governed by the Hamiltonian ^ 

H2DIM = «(1 + + 2a cos ki — 2 cos k2)a^ + 2a(sin kiai + sin /C2C2) 

with a = -\/2 — 1- This Hamiltonian is similar to Hi as it has only one node at /c = 0. 



IV. CONCLUSIONS 

The previous discussion is based on the one-particle Green's function without 
referring to an effective field theory. This approach is motivated by the aim to avoid 
additional symmetries which appear by the introduction of new fields and which 
are not related to the original Hamiltonian or the one-particle Green's function. In 
the case of the Hamiltonian H, regardless of the discretization of Vj, there are two 
discrete symmetries: 
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H ^ -a-sHcxs (4) 

for the ensemble and, therefore, for the average one-particle Green's function. More- 
over, there is a symmetry under the discrete transformation 

H -aiH^au (5) 



which defines class D of Ref. [|T9[. It holds for each realization of the model with 
Hamiltonian H. hnnn and /i^„„, both break symmetry whereas only /?.„„„ breaks 
symmetry (P). 

From the symmetry point of view alone it is not entirely clear under which 
conditions a vanishing average DOS at ii^ = exist. The competition of different 
nodes (i.e. different length scales) cannot be described only by global symmetries 
but requires more detailed knowledge. However, there is a simple argument in terms 
of the hopping expansion which indicates that the number of independent random 
terms in the Dirac Hamiltonian H is crucial for the behavior of the average DOS 
around E = 0. Taking the zero-dimensional limit of H, i.e. the leading order of the 
hopping expansion, there is a power law 

{p{E)) Po\Er, 

where a = (random mass) and a = 1 (two-component random vector potential). 
The latter, of course, violates also the symmetry (||). 

The non-vanishing DOS for the class D model in d = 2 with a diffusive behavior 
was recently discussed by Senthil and Fisher and Read and Green in terms 



of the RG flow of a non-linear sigma model. Although this approach was criticized 
by Bocquet et al. [T^, its result is in agreement with that obtained with a different 



approach . The appearence of the intermediate phase is quite natural for the 2D 
Dirac fermions, since the pure model has a singular "metallic" phase a.t E = m = 
with conductivity /hn. This phase is robust against a random vector potential, 
where the value of the conductivity remains unchanged [Q. A random mass has 
apparently a stronger effect because it reduces the conductivity at m = by a 
factor 1/(1 -|- (yf/vr), where g is the variance of the random mass [Q. Moreover, it 
broadens the singular phase at m = to an interval — < m < rric with non- 
vanishing DOS. In terms of the random bond Ising model the non-vanishing DOS 
reflects the existence of the Griffiths-McCoy- Wu phase, which cannot be seen in the 
perturbative RG approach 



The occurence of vortex-like excitations might be an important effect, as sug- 
gested in Refs. [0,|^,|1^. Using the vr flux discretization, these vortices can be 



created by local edge currents in areas where the sign of the random mass changes: 
An area with a positive mass has a positive Hall conductivity whereas the sur- 
rounding area with negative mass has a Hall conductivity with opposite sign. The 
resulting edge currents can have a long-range interaction. An effective model for this 
behavior can be found in terms of the Q matrix fleld theory of Ref. |jl2|], in which 



the random mass is replaced by a matrix fleld. Details will be published separately. 

In conclusion, the change of the discretization has a strong effect on the node 
structure of the Dirac Hamiltonian. The correct discretization is determined by 
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the effective gauge field which the Dirac fermions experience. In the case of Dirac 
fermions with random mass, however, the average DOS at low energies is relatively 
robust against the change of the discretization. In particular, the k ^ modes 
have a substantial contribution to the average DOS. This indicates that also the 
intermediate metallic phase of the Dirac fermions with random mass and vr flux 
discretization should be robust under a change of the discretization. 
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FIG. 1. Simple plaquettes are created from the Hamiltonian Hq with nearest neighbor 
discretization. They enclose a flux (f) = tt which gives e**^ = —1. The sublattice gauge 
transformation H — > HSa^ does not change this property. 
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